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ABSTRACT 

We present Hubble Space Telescope NICMOS i/-band imaging of 33 ultraluminous infrared 
galaxies (ULIRGs) at z ~2 that were selected from the 24 fim catalog of the Spitzer Extragalactic 
First Look Survey. The images reveal that at least 17 of the 33 objects are associated with 
interactions. Up to one-fifth of the sources in our sample could be minor mergers, whereas 
only two systems are merging binaries with luminosity ratio <3:1, which is characteristic of 
local ULIRGs. The rest-frame optical luminosities of the sources are of the order 10 10 -10 11 
Lq and their effective radii range from 1.4 to 4.9 kpc. The most compact sources are either 
those with a strong active nucleus continuum or those with a heavy obscuration in the mid- 
infrared regime, as determined from Spitzer Infrared Spectrograph data. The luminosity of 
the 7.7/im feature produced by polycyclic aromatic hydrocarbon molecules varies significantly 
among compact systems, whereas it is typically large for extended systems. A bulge-to-disk 
decomposition performed for the six brightest (r?iif<20) sources in our sample indicates that 
they are best fit by disk-like profiles with small or negligible bulges, unlike the bulge-dominated 
remnants of local ULIRGs. Our results provide evidence that the interactions associated with 
ultraluminous infrared activity at z ~2 can differ from those at z ~0. 
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The excess of infrared (IR) emission in ex- 
tragalactic sources is widely believed to be trig- 
gered most efficiently b y gas-rich galactic merg- 
ers or interact i ons (e .g, lLarson &: Tinslev" 1978t 
Sanders et al. | Il988al) . Dynamical instabilities 
are required to drive large concentrations of 
gas at small radii and to lead to a high rate 
of star-formation and an excess of IR emis- 
sion (> 10 11 Lq). The number density of 
both IR-bright sources and interactions evolves 
strongly with redshift z, increasing to a redshift 
of at le ast 1 (lElbaz et al. I [2003 iLe Floc'h et al ' 



Observations of ultraluminous IR galaxies 
(ULIRGs; L 8 _i oo M m > W 12 L Q ) in the lo- 
cal Universe indicate that ULIRGs represent a 
well-defined merger category. They are trig- 
gered by the merger of two gas-rich galax- 
ies whose typical gas frac tions are ~10% of 
their dynamical masses (jDownes fc Solomon 
19981 ). The merging components have compa- 
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Perez-Gonzalez et al 
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2006; iKartaltepe et al. 



20051 
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rabl e masses flDasvra et al. I l2006af) or lum inosi- 
ties (|Kim et al. 1120021 : iMurphv et al. Ill996l) . The 
fairly specific origin of ULIRGs is thought to 
be related to their high phase-space gas densi- 
ties, which are achieved through major merg- 
ers (e.g., iBarnes I Il992h . The merger typi- 



However, it is still 
unclear whether the interaction mechanisms that 
lead to high IR outputs are the same in the local 
Universe and at high redshifts. 



call y leads to the formatio n of elliptical galax- 
ies (jjosephfe Wright Ill985h . The remnants often 
have light profiles that follow an r 1//4 distribu- 



tion (jKormendv fc Sanders Ill992t IVeilleux et al 
2002), fall on the velocity dispersion — cffcc- 
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tive radius iL ff proj ec tion dGenzel et al. _ 2001 



, . 'p i 

Tacconi et all 120021 iDasvra et al. I l2006bh "of 



the fundamental plane of early-type galaxie s 
dDressler et al. Ill987t [Diorgovski fc Davis Ill987l). 



and h ave a bulge-to-disk ratio of ~ 1.5 (jVeilleux et al. 
20061) . 

The trigger mechanisms and end products 
of lower luminosity systems, such as the lo- 
cal luminous IR galaxies (LIRGs; Ls-iooo/jm ~ 
10 11 — 10 12 L Q ), are not constrained as well as 
those of local ULIRG s . The imaging analysis of 
Alonso-Herrero et al. I ( 20061 ) indicates that most 
local LIRGs have prominent spiral patterns, and 
of those, a nonnegligible number are weakly in- 
teracting or even isolated systems. Ilshida I (|2004l ) 
found that several minor mergers or fly-by inter- 
actions exist at the low-luminosity end of local 
LIRGs. However, as the luminosity increases 
so does the fraction of major mergers. Some 
LIRGs need to have the same origin as ULIRGs 
since every ULIRG has to pass through a LIRG 
phase before and afte r its main starburst episode 



(jMurphv et al. Il200ll ) 



At high redshifts the situation is not as well 
defined, mostly because the merging galaxies are 
on average more gas rich than those in the lo- 
cal Universe. The amount of gas in, e.g., sub- 
millimeter galaxies i s ~40% of their total mass 
( Tacconi et al. If2006h . It has been suggested that 
a variety of galaxy types can be produced by 
mergers of high gas fr actio ns. The simulations of 
Springel et al. I ( 20051 ) and Springel fc Hernquist I 
(|2005 ) indicate that the merger outcome can vary 
between a red and dead elliptical or a spiral galaxy, 
depending on how much gas is converted into 
stars o r expelled from the galaxy. Hammer et al. I 
(|2005l ) suggest that most local spiral galaxies may 
have undergone mergers within the last 8 Gyr by 
invoking star-formation duration and mass assem- 
bly arguments. 

The goal of this paper is to investigate what 
processes lead to an ultraluminous IR activity at 
redshifts 1.5< z <2.5. Could high-z ULIRGs be 
triggered in an analogous manner to local ULIRGs 
or to local LIRGs? If triggered by interactions, 
what are their progenitors and their end prod- 
ucts? To address these questions, we requested 
observations of a sample of 35 ULIRGs at z ~2 
using the Near Infrared Camera and Multi-Object 
Spectrometer (NICMOS) instrument on board the 



Hubble Space Telescope (HST) . 

The rest of this paper is organized as follows. 
In §[2] and §[3] we present the sample and the data 
reduction used for the analysis in § SJ The presen- 
tation of the results in §[5] is followed by a compar- 
ison between z^2 and local ULIRGs in §[6] and the 
conclusions in § [71 Throughout this paper we use 
a ACDM cosmology with H o =70 km s" 1 Mpc" 1 , 
Q ro =0.3, and fi A =6.7. 

2. The Sample 

Our sources are drawn from a sample of IR- 
bright galaxies detected in the 24 jum mosaic of 
the 3.7 deg 2 Spitzer Extragalacti c First Look Sur- 
vey (XFLS; iFadda et al.Tl2006j) . The selected 



sources have 24 u m fluxes greater than 0.9 mJy 
( Yan et al~l I2005T I . Two extra criteria were im- 



posed to ensure that the sample includes high-z 
starburst galaxies: the flux ratios log(^/ l/ [24 /im]/ 
vf v [8 H) and logO/424 H/ vf v [6440 A]) were 
required to be greater than or equal to 0.5 and 
1, respectively. The use of these threshold val- 
ues is based on template spectr al energy distribu - 
tions (SEDs) of local starbursts (|Yan et al. II2004I ). 
The high IR-to-optical flux ratio preferentially se- 
lects starbursts at z > 0.6, since the observed R- 
band magnitudes sample the SEDs shortward of 
the 4000 A break. The high 24/8 fim flux ratio 
can either pick strong starbursts at 1.5 < z < 3 
whose 6.2, 7.7, and 8.6 /jm polycyclic aromatic 
hydrocarbon (PAH) features are shifted in the 24 
/j,m beam or obscured, reddened active galactic 
nucleus (AGN) con tinua (Yan et al. 2004, 2007; 
Brandl et af~ll2006h . In total, 52 XFLS sources 
that satisfy these criteria wer e observed with the 
Infrared Spectrograph (IRS; iHouck et al~l 12004 ) 
on board Spitzer. The mid-infrared (MIR) spec- 
tra yielde d redshifts for 4 7 sources in the range 
[0.61,3.2] (lYan et al. Il2007h . For the NICMOS ob- 
servations, we chose to observe only the sources at 
z > 1.5 so that we can perform an imaging anal- 
ysis at comparable rest-frame wavelengths for all 
sources. Of the 35 sources at z > 1.5, 33 were 
successfully observed. These 33 sources are pre- 
sented in Table [T] and their redshift distribution is 
presented in Figure [TJ 

From the IR point of view, most of the sources 
in this sample arc ULIRGs. The rest-frame 14 
/im luminosity, L14, which is based on spectral 
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fitting of IRS data by ISaiina et al. ( 2007 ). is 
greater than 10 12 L Q for 27 sources, indicating 
that their bolometric IR luminosities, Lir, are 
in the range of ULIRGs or higher. The bolo- 
metric IR luminosities are computed, when pos- 
sible, by fitting the far-infrared (FIR) SEDs of 
the sources using Spitzer MIPS 70 and 160 fim 
fluxes, and Max-Planck Millimeter Bolometer 1.2 
mm fluxes (Sajina et al. 2008; see Table []}. 
The SED fitting indicates that three more sources 
have Lir values that exceed the ULIRG threshold. 
On average, the IR luminosity of these sources is 
subst antially greater th an that of local ULIRGs 
(e.g., iKim et al. I l2002h and is close to the 10 13 
Lq limit, which defines the hyper luminous IR 
galaxies. The remaining three sources, namely 
MIPS 8493, MIPS 16133, and MIPS 22661 are also 
likely to be ULIRGs given their high L14 values 
and Lir, upper limits (Sajina et al. 2008). 

The 7-38 /im IRS spectra revealed that the 
IR emission arises to comparable levels from in- 
tense star formati on, AGN activity , or a combi- 
nation of the two ( Yan et al. 1120071 ). About one- 



third of the sample has a strong silicate absorp- 
tion feature at 9.7 /im, which implies that the 
main emission mechanism is compact and heav- 
ily obscured; it c an either be a nucl ear starburst 
and/or an AGN (| Spoon et al. Il2007h . Therefore, 



the IR emission of high-z ULIRGs has a multiple 
origin, similar to that of local ULIR Gs as revealed 



by Infrared Space Observatory (ISO; lGenzel et al. 



1998; iRjgopoulou et al. I fl9 991 and Spitzer spec- 
troscopic data ([Armus et al. 1120071 ). 

3. Observations and data reduction 

The observations of this program (ID 10858) 
were carried out using the NIC2 camera of the 
NICMOS instrument. The pixel scale of NIC2 cor- 
responds to 0'.'076 and 0'.'075 along the x and y 
axis respectively. The resolution of the images is 
high, 0V158, as measured from the full width half 
maximum (FWHM) of stellar light profiles. The 
images were taken in the i?-band using the F160W 
filter, which is centered at a wavelength of 1.6 /im. 
The wavelength that corresponds to the center of 
the filter at the rest frame of each source is in the 
range [4310, 6150] A. The width of the F160W 
filter at the rest frame of each source is presented 
in Table Q] The detector read-out mode was set 



to "multiaccum" with 23 or 24 consecutive, non- 
destructive reads leading to an individual expo- 
sure. The observations were performed using a spi- 
ral four-point dither pattern per orbit to optimize 
the removal of bad pixels. The total on-source 
integration time for the four dither positions was 
either 2688 or 2560 s. The number of orbits per 
source was determined by the source's Spitzer In- 
frared Array Camera (IRAC; iFazio et al~l 120041) 
3.6 fim flux. The sources with flux density greater 
than the detection lim it of the XFLS catalog (20 
/itJv: lLacv et al. 2005) were typically observed for 
a single orbit. The sources that are fainter than 
this threshold were observed for two orbits. The 
total on-source integration times are presented in 
Table [TJ 

The reduction and calibration of the NICMOS 
images were performed using standard HST IRAF 
and IDL routines and foll owing the prescri ptions 
of the NICMOS manual ([Barker et al. 1120061) . The 
pipeline calnica function was used to remove the 
zeroth detector read-out image from the non- 
destructive consecutive read-out images, and to 
correct for detector nonlinearities, dark currents, 
and bar effects. The same task was also used to 
flat field and combine the consecutive read-out im- 
ages into a single image, computing its statisti- 
cal error. We then used the IRAF task pedsky 
to bring the quadrants of the detector to similar 
bias levels and to apply a zeroth order removal 
of the sky background. Following pedsky, we cor- 
rected for South Atlantic anomaly effects using the 
IDL procedure saa clean. To correct for system- 
atic sky residuals and enhance the signal-to-noise 
ratio of our background-limited observations, we 
constructed a residual sky image by median com- 
bining the images at the different dither positions. 
We then removed the sky residuals from each of 
these images, shifted them into their common ref- 
erence position, and removed cosmic rays and bad 
pixels. The fluxes of the bad pixels were replaced 
with their average value in all other useful dither 
positions. To enhance the resolution of our im- 
ages, we resampled the pixels to half their origi- 
nal size using the IRAF routine drizzle. This rou- 
tine simultaneously corrects for detector geomet- 
ric distortions. The resulting images were weight- 
averaged using their uncertainty images, which 
were also corrected for shifting and drizzling trans- 
formations. The weight-averaged image and its fi- 
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nal uncertainty image were drizzled back to their 
initial pixel scale. The final images are presented 
in Figure [21 

4. Data analysis 

Our images show that a significant fraction of 
sources have irregular patterns indicative of in- 
teractions. For the proper photometric calibra- 
tion of such sources, we need to ensure that all 
of their light will be accounted for, even if it is 
distributed in asymmetric structures (e.g., tidal 
tails). For this reason, we performed a basic pho- 
tometric analysis by fitting a series of isophotal 
ellips es to each galaxy with t he sextractor soft- 



ware ( Bertin Sz Arnouts 19961 ). Using sextractor, 



we first ensured the detection of each source at 
a 3cr level and then integrated its flux down to 
la level. Sextractor measured the flux, the ef- 
fective (half-light) radius, the ellipticity, and the 
position angle of each galaxy within that detec- 
tion limit. We converted the integrated galaxy 
counts per second to an apparent Vega magnitude 
using the NICMOS da ta calibration prescription 
(|Dickinson et al. 1120021 ) 



m = -2.51og(/ CPS x PHOTFNU/ /„ [Vega]), 

(1) 

where feps is the source flux in counts per second, 
PHOTFNU is the conversion factor from counts 
per second to janskys, and /^(Vega) is the Vega 
flux in janskys for the given NICMOS camera and 
filter passband. PHOTFNU is computed by cal- 
nica and equals 1.49585 xlO -6 Jy s count -1 . For 
NIC2 at F160W, /„(Vega) is 1043.6 Jy. 

Following the global photometric characteriza- 
tion of all the sources in our sample, we decom- 
posed the brightest syste ms into a bulge an d a disk 
component using galfit ( Peng et al. 1120021 ). Galfit 
implements a \ 2 minimization to decompose the 
light profile of a source into a (set of) predefined 
function(s). We selected the Sersic profile, whose 
intensity scales with radius as exp([r/i? ff] 1 /™). 
We used two Sersic components, one with n=4 
and another with n=l, corresponding to a de Vau- 
couleurs and an exponential disk profile respec- 
tively. Galfit convolves both components with the 



instrumental point-spread function (PSF) before 
simultaneously fitting them to the observations. 
We also performed a PSF versus free Sersic in- 
dex decomposition as a consistency check. We 
simulated the required NIC2 P SF with the Tiny- 
Tim code (|Krist fc Hook Ill997| ) using a pixel scale 
identical to that of the observations. To perform 
the x 2 minimization, galfit requires an error image. 
For this purpose, we provided the final uncertainty 
image. Galfit also requires initial guesses for the 
parameters under examination; we used the out- 
put parameters of sextractor as input values for 
the magnitude, the half-light radius, the position 
angle, and the ellipticity of each source. Finally, 
we masked out all other bright sources within each 
image so that galfit can properly compute its back- 
ground level. 

5. Results 

5.1. Rest-frame Optical Galaxy Morpholo- 
gies And Luminosities 

The mean and median i7-band (or rest-frame 
optical) magnitudes of the sources in our sample 
are m#=20.89 and mjy=21.07 with a dispersion 
of 0.98 in Vega magnitudes (Tabled Fig. ED . For 
comparison, one of the brightest ULIRGs in the 
local Un iverse, Mrk 1014, would have m#=20.7 
at z=2 dKim et al. 1120021 ). The distribution of 



2 Galfit is thought of as the optimal code for the extraction 
of morphologica l parameters of unrelaxed galaxies at high 



redshifts (Haussler et al. 2007). 



the half-light radii is also presented in Figure [3] 
The mean and median values of the R c g distribu- 
tion equal 2.66 and 2.43 kpc, respectively, with a 
standard deviation of 0.80 kpc. The R — H and 
H— IRAC(3.6 fim) color distributions (Fig. [3]) have 
mean values of 1.86 and 2.20 and median values of 
1.64 and 2.38, respectively. Their standard devia- 
tions equal 1.06 and 0.74. The surface-brightness 
Hh distribution has a mean, median, and standard 
deviation of 20.33, 20.62, and 1.06 mag arcsec~ 2 
(Table [2]). We computed the optical luminosity 
-^opt of each source at the rest-frame wavelength 
that corresponds to the center of the F160W fil- 
ter (1.6 /im) by multiplying its flux with the fre- 
quency at the center of the filter. We find that the 
rest-frame optical luminosities range between 1.70 
xlO 10 L and 9.46 xlO 11 L Q (Table [2]). The av- 
erage and median luminosity values are 1.12 x 10 11 
Lq and 7.83 xlO 10 Lq with a wide dispersion 
that is equal to 0.37 dex when measured on a log- 
arithmic scale. The rest-frame L14/ L opt ratios 
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are in the range [4,218]. For the six sources with 
determined IR luminosities, the LmJ L op t ratios 
range between 50 and 281. 

The NIC2 images show that 2 of the 33 sources 
(MIPS 16059 and MIPS 22530) have two merging 
components of luminosity ratio <3:1. There are 
three more systems, MIPS 289, MIPS 8327, and 
MIPS 15928, with a secondary faint component 
inside the envelope of their primary nucleus. The 
faint component could either be another nucleus 
or a star-forming cluster. Other signs of interac- 
tions, such as perturbed morphologies, stellar fans 
(MIPS 22558), and tidal tails, exist in 12 more 
sources. MIPS 8196 has a bright tail of 19.3±3.7 
kpc. In total, 52% of the sources in our sample 
appear to be involved in interactions. 

This fraction is likely to be a lower limit, given 
observational limitations. Wide binaries that lack 
prominent signs of interaction or close binaries 
that are unresolved may exist in our sample. In 
the case of unresolved systems, their maximum 
projected nuclear separation d n is determined 
by the instrumental angular resolution, which 
equals 2.1 pixels or 1.34 and 1.28 kpc at redshifts 
of 1.5 and 2.5, respectively. Such a resolution 
is comparable to that of ground-based observa- 
tions that were used to study large numbers of 



ULIRGs in the local Univ erse (e.g., iKim et al. 



2002; IVeilleux et al. 1 12002). which was ~1.8 kpc 



at z=0.1. Higher angular resolution imaging stud- 
ies have shown that a small number of local 
ULIRGs are likely to have two nuclei separated 
by less than ~1 5 kpc (Scoville et al 1998 ; 2000; 
Bushouse etai~ll2002t iDasvra et al. ll2006tJ ). Flux 
or surface-brightness detection limitations could 
also contribute to this effect. For a single orbit, 
the flux 3<7 detection limit of our images corre- 
sponds to m#=26.07 for a point source and to 
m#=24.51 for an extended source with a radius 
equal to the average radius of the sources in our 
sample. This value implies that for a primary nu- 
cleus of 771^=22.00, the secondary nucleus needs 
to have a Am <2.5 in order to be detected. Such 
a magnitude difference corresponds to a lumi- 
nosity ratio of 10:1. Therefore, small compan- 
ions in mergers of high luminosity ratio could 
be undetected. Moreover, low surface brightness 
structures that have been diluted into the sky 
background could exist. The surface brightness 
3(J detection limit of our observations is 22.59 



mag arcsec -2 for a single orbit. Tidal tails in 
local ULIRGs have B-band surface brightness val- 
ues between 21.5 mag arcsec -2 , e.g., in IRAS 
22491-1808 and IRAS 14348-1 447, and 23.5 mag 
arcsec" 2 , e.g., in UGC 5 101 (jSurace fc Sanders 
120001 ISurace et al. Il200(l . 



To further investigate the effects of surface 
brightness dimming on our results, we created a 
model that simulates z=2 galaxies from local tem- 
plates. The chosen templates, the methodology, 
and a detailed description of the results can be 
found in the Appendix. We find that the initial 
light profile is usually maintained for sources with 
777/f=19. Extended structures are typically diluted 
into the continuum for sources with mjy>20. It 
is therefore likely that a nonnegligible fraction of 
our sources could have diluted low surface bright- 
ness features. Up to half of the flux can be di- 
luted into noise for sources with mu—21, leading 
to effective radii that can be underestimated by 
up to a factor of 2. Since the average and me- 
dian effective radius of the local ULIRGs in the 
1 Jy sample are 4.85 and 4.70 kpc, respectively 
(Vcil leux et al. 2002). he intrinsic radial extents 
of these z^2 ULIRGs are likely to be comparable 
to those of local ULIRGs. 

In all computations we have ignored any possi- 
ble extinction or fc-corrections in t he rest-frame 
wavele ngth range [4200,6200] A. iTavlor et a~ 



(|2005l ). who performed imaging of a sample of 
142 late-type and irregular or interacting galax- 
ies, found that their average B — V color is 0.5, 
corresponding to a 37% difference in the flux of 
the two bands. This difference in flux does not 
necessarily correspond to a systematic difference 
in R e g from the B to the V band. Often, large 
color gradients between t he center and the o uter 
regions of ULIRGs exist (jSurace et al. 1119981 ) be- 
cause the geometry of the o bscuring medium is 



complex or highly nucleated (jSoifer et al. 1 12002) 



due to dynamical instabilities. Moreover, account- 
ing and correcting for B — V color differences is 
not straightforward since the correction can sig- 
nificantly deviate from its average value in indi- 
vidual objects. Spectroscopy of local ULIRGs has 
shown that the continuum slope in the optical 
wavelengths varies from source to s ource because 



of diff erent AGN-emission strength (jArmus et al 



19891 Veilleux et al. 1995; 1999). The broad-band 



colors of nearby filters may also significantly dif- 
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fer due to bright emission lines within the filter 
bandpasses. 

5.2. Bulge-to-Disk Decomposition Of Bright 
Sources 

We used galfit to find the stellar distribution in 
z ~ 2 ULIRGs. We performed this analysis only 
for the six systems with ijih< 20, because the light 
profiles of sources with mn >20 can be diluted in 
the background, leading to unreliable decomposi- 
tion results. To infer this reliability limit, we ran 
galfit for two local disk- and bulge- dominated sys- 
tems and their simulated z=2 images. This anal- 
ysis is presented in detail in the Appendix. 

The bulge-to-disk decomposition results are 
presented in Figure |4] and in Table [3] The resid- 
ual flux of the decomposition, which on average 
corresponds to 4.44% of the initial galaxy flux, 
is also presented in Table [3] The decomposition 
reveals that if bulges indeed exist in the sources 
that we analyzed, they are small. Among the six 
sources, four have a possibly unresolved bulge. 
This means either that the bulge is real but its 
size is smaller than our resolution element or that 
it does not exist and galfit identifies a compact 
nuclear starburst or an AGN as a bulge. The 
mean bulge-to-disk ratio of these sources is 0.62. 
The median is 0.55, and it is significantly smaller 
than that of local single-nu cleus ULIRGs, which is 



1.51 (IVeilleux et al. 1120061) . The measured bulge- 



to-disk ratios may be upper limits of their in- 
trinsic values, since the surface brightness dim- 
ming typically affects more the disk- than the 
bulge- dominated light profile distributions (see 
the Appendix). These results imply that most 
of these systems had not formed big bulges as 
of -10 Gyr ago. The exception is MIPS 8196, 
which is the brightest system in our sample with 
vL v — 9.46 x 10 11 Lq and a bulge-to-disk ratio 
>1. 

Sin ce a small bulge can be degenerate with a 
PSF (S immons fc Urrv 1 120081 ) in decomposition 
algorithms, we also performed a PSF versus free 
Sersic index n decomposition to check the con- 
sistency of our results. We find that the effective 
radii of the free-index component agree on average 
to 8.53% with those computed using SExtractor. 
We also find that the PSF flux is always fainter 
than the galaxy flux by more than a factor of 3 
(or 1.22 mag), and that in five of the six cases the 



index n is close to unity. This result confirms that 
the systems are disk-dominatecH. However, this 
analysis needs to be performed on large samples 
to investigate how frequent or efficient the forma- 
tion of bulges by ULIRGs is at z=2. 

A simila r bulge-to-dis k deco mposition was per- 
formed by IZheng et al~l (|2004h on lower redshift 
(z <1.2) LIRGs and ULIRGs. Their sources were 
detected at 15 /im using ISO Ca mera data in 



the C anada-France Redshift Survey (jFlores et al 
1999). Their sources were observed with the HST 



Wide Fiel d Planetary Camera 2 in the B,V , and 
I bands. IZheng et al. I (J200J) also find that the 
six single-nucleus ULIRGs in th eir sample likewise 
have s mall or negligible bulges. iMelbourne et al.l 
(|2005l ) performed a morphological analysis of IR- 
bright sources up to z—1 in the Great Observato- 
ries Origins Deep Survey (GOODS) 24 /im image 
(R. R. Chary et al. 2008, in preparation) and found 
that z—l LIRGs are mostly spirals or peculiar. In 
the combined GOODS and Galaxy Evo lution from 



Morph ology and SED (GEMS) fields, iBell et al. 



(|2005l ) find that approximately 50% of the sources 
are spiral galaxies at z=0.7. 

5.3. Rest-frame Optical vs MIR Galaxy 
Properties 

In Figure [5] we examine the dependence of the 
rest-frame optical morphological parameters of the 
sources in our sample on their MIR spectral classi- 
fication. The optical magnitudes and luminosities 
are plotted as a function of the effective radii for 
the sources with well-defined MIR spectral types 
(see Table [1]). Both the magnitudes and the lumi- 
nosities are statistically indistinguishable for all 
types of sources. However, the various systems 
segregate along the effective radius axis. AGN- 
dominated sources are compact systems with R c g 
of —2 kpc, whereas PAH-dominated, starburst 
sources are as extended as 5 kpc. The obscured 
sources have a large range of effective radii, com- 
prising both compact and extended systems. 

In Figure [6] we plot the 7.7 /xm PAH feature 
luminosity, I7 7, a nd equivalent width, EW7.7 
(jSaiina et al. 1 120071 ). as a function of the rest- 



3 We note that disk-like systems at z ~2 should not neces- 
sarily be thought of as spiral galaxies. Disks mostly refer 
to all smooth or somewhat clumpy light distributions that 
follow an exponential profile. 
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frame optical effective radius. In this plot, we do 
not distinguish between sources of different MIR 
spectral types. We find that high PAH luminosi- 
ties appear in both compact and extended sources 
but the lowest values of L7.7 appear in compact 
systems. The 7.7 [im equivalent width does not 
correlate with the optical extent of the system, 
but the sources with the highest EWs are ex- 
tended. For compact objects, this result proba- 
bly reflects the dilution of the PAH features by 
the AGN continuum, since compact sources have 
on average a stronger AGN component in both IR 
and opti cal wavelengths than most extended sys - 



Spectrograph (STIS) ob servations of 1< z <3 



li- 
tems do ([Tran et al. 1120011 IVeilleux et al. 1 120021) . 



This result could be a confirmation of a scenario 
in which ULIRGs with warm MIR/FIR SEDs are 
more AGN-domin ated than ULIRGs w ith cold 
MIR/FIR SEDs (ISanders et al. I Il988br ). Still, 
compact sources with high PAH luminosity exist. 

In Figure we plot the MIR optical depth, 
t, as a function of the rest-frame optical effective 
radius and luminosity of each source. The value 
of t is der ived from the SiO absorption feature 
at 9.7 /xm ([Sauna et al. I l2007h . The spatial dis- 
tribution of the silicates that are responsible for 
the MIR obscuration is not correlated with the 
rest-frame optical stellar light profile. The optical 
depth of the MIR obscuring medium does not cor- 
relate with the stellar luminosity either: obscured 
systems in the MIR can be bright in the optical 
wavelengths (e.g., MIPS 8196). This result could 
imply t hat the MIR obscur i ng medium is circu m- 



nuclear ([Spoon et al. 1120071 : ISaiina et al. 1120071 ). 



6. Discussion: Trigger mechanisms of z~2 
vs local ULIRGs 

A qualitative morphological classification of the 
ULIRGs in our sample indicates that there is a 
significant fraction of sources that appear to be 
quiescent or compact. The number of sources 
that have no prominent signs of interaction is 
48%. This fraction is high in comparison with 
the local ULIRG samples, in which only 2% of 



the sources are not distorted ([Kim et al. 1 120021 : 
Veilleux et al. 20061 ). The high fraction of qui- 



escent sources could be attributed to flux or sur- 
face brightness limitations of our observations (see 
However, HST NICMOS, Advanced Camera 
for Surveys (ACS), and Space Telescope Imaging 



submillimeter galaxies ( Chapman et al. I 120041 



Swinbank et al. 



20061 ) indicate that ~80% of 
these sources have signs of interaction. Although 
color effects may change this fraction, it is possible 
that the differences could alternatively be specific 
to the selection criteria of each sample. Since our 
objects are brighter than 0.9 mJy at 24 /im, they 
have warmer MIR/FIR SEDs than submillime- 
ter galaxies (Sajina et al. 2008). According to a 
scenario in which a warm ULIRG phase follows 
a cold ULIRG phase in the sa me evolutionary se- 
quence ([Sanders et al. 1l988bl ). our sources appear 
to be less perturbed than submillimeter galaxies 
because they could be closer to dynamical relax- 
ation. In this picture, the submillimeter galaxies 
could be the cold, binary mergers and the 24 /xm- 
selected ULIRGs could be their warm, merged 
counterpar ts. This possibility i s supported by the 



findings of iBridge et al. I (|2007[ ). who showed that 



the likelihood of a merger to be in a close pair is 
5 times higher for sources brighter than 0.1 mJy 
at 24 /jm than for sources below this limit. 

To investigate for (dis)similarities in the prop- 
erties of z^2 and local ULIRGs, we compared the 
number of on-going major mergert@ in our sam- 
ple and i n the entire 1 Jy catalog of mainly local 
ULIRGs dKim et al. 1120021 ). Of the 118 sources in 
the I Jy catalog, 40 (or 34%) are binaries with an 
optical (B— to R— band) lumin osity ratio <3:1, or 



Am<1.2, and 1.5< d n <30 kp<H (|Kim et al. II2002I: 
Surace et al~ll2000l ). In the z~2 sample, 2 (or 6%) 



of the 33 systems (MIPS 16059 and MIPS 22530) 
are currently interacting binaries with Ato<1.2 
and 1.5<d n <30 kpc. We ran a Monte Carlo code 
to quantify the significance in the difference of 
these fractions. We created a parent population of 
118000 sources, 40000 of which are binaries that 
follow the d n and Am distribution of the I Jy sam- 



4 We note that all comparisons between the local and the 
z~2 ULIRGs presented in this Section are based on the rest- 
frame optical luminosity ratios of the interacting galaxies. 
Since this ratio could be affected by extinction, it could 
deviate from its intrinsic value. Local ULIRGs that appear 
to be minor mergers in optical images but major mergers 
in near-IR images and vice versa exist l|Kim et al. 2002: 
ISurace et al~l 120001) . The strength of this effect can vary 
significantly from source to source, depending on the geo- 
metric distribution of its obscuring medium. 

5 We chose our fiducial nuclear separation threshold to be 
30 kpc since there is only one binary system in the 1 Jy 
sample with d n that exceeds this value. 
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pie binaries. We randomly chose 33 sources and 
computed the probability p that a number ns of 
the 33 sources meet the criteria 1.5<d n <30 kpc 
and Am<1.2. We iterated the procedure for 10 5 
times and found that p is insignificant (0.014%) for 
n s =2. Therefore, our sample statistically differs 
from the entire 1 Jy sample. Four more sources in 
our sample (MIPS 464, MIPS 15840, MIPS 15880, 
and MIPS 22651) have neighbours that satisfy our 
nuclear separation and luminosity ratio criteria. 
However, they have no signs of interaction with 
their neighbours. If we assumed all these sources 
to be interacting, the maximum fraction of major- 
merger pairs in our sample would increase to 18%. 
In that case, the probability p would be 3.0%, 
which is significant but still small. 

We then examined whether the sources in 
our sample are the analogs of the warmest lo- 
cal ULIRGs by comparing the number of on-going 
major mergers in both our sample and in a local 
luminosity-matched control sample. We selected 
the control sample as the 1 Jy catalog sou rces 
that have L m > 10 125 L Q (|Kim et al. ll2002l ) for 
the cosmological parameters adopted in this pa- 
per. The local control sample has 16 sources that 



gered by such mergers exist, e.g., IRAS 12127- 



1412 (|Kim et al. I l2002h . Up to one-fifth of the 



typically have warm MIR SEDs (e.g.. iTran et al 



2001; lArmus et al. 1120071) . Two or 13% of these 



sources are binaries with Ato<1.2 and 1.5<e?„<30 
kpc. This fraction of binary systems is statisti- 
cally indistinguishable from that in our sample 
(6% -18%). Therefore, the z~2 ULIRGs mor- 
phologically resemble their warmest local analogs. 
However, they differ from their local analogs in 
their MIR spectral properties. IRS spectroscopy 
indicates that the majority of the most IR lumi- 
nous local sources have a strong AGN component 
in the MIR (L. Armus et al. 2008, in prepara- 
tion), unlike the z^2 sample in which systems 
with strong starburst-related features exist. One 
such example is MIPS 22530, which has a MIR 
spectrum with strong PAH emission and Lir = 
8 x 10 12 Lq. We conclude that there is no lo- 
cal ULIRG (sub)sample that is likely to comprise 
sources with simultaneous compact morphologies, 
high IR luminosities, and strong starburst-related 
features, which are characteristic of these z~2 
ULIRGs. 

In the local Universe, minor mergers are asso- 
ciated with LIRGs instead of ULIRGs, although 
few atypical examples of ULIRGs that are trig- 



ULIRGs in our sample could be associated with 
on-going minor mergers. Six sources (MIPS 78, 
MIPS 180, MIPS 8242, MIPS 8342, MIPS 8493, 
and MIPS 16113) have low surface brightness 
neighbours within a radius of 10 kpc and could be 
close interacting binaries. We used galaxy number 
counts N g per square degree and magnitude to find 
the probability of randomly having a faint back- 
ground or foreground source within a projected 
distance of 10 kpc or 1'.'2 at 2=2. We constructed 
a Monte Carlo simulation to randomly distribute 
N g galaxies over 1 deg 2 and to compute the pro- 
jected separation of each source from its nearest 
neighbour. For ^=150000 galaxies in the range 
m H =22.5 to 23.5 (Metcalfe et al. II2006I ). we found 
that 5 out of every 100 sources have a neighbour 
within the given separation due to projection ef- 
fects. We then selected 33 sources from this parent 
population and computed the probability that at 
least six galaxies have neighbours due to projec- 
tion effects. We iterated the procedure for 10 5 
times and found that the probability is 0.52%, 
which indicates that some of these sources are 
likely involved in on-going minor mergers. 

In £j5j we showed that five of the six brightest 
sources in our sample have a small Sersic index 
(< 1.35) and bulge-to-disk ratio (< 0.57). Such 
systems can evolve to disk galaxies with small 
bulges, unlike local ULIRGs. Regardless of how 
representative this result is for the majority of z=2 
ULIRGs, it provides further evidence that the dy- 
namical configuration of IR-bright sources can dif- 
fer at z~2 and in the local Universe. 

A plausible explanation regarding why some 
sources differ from their local analogs is that the 
z^2 systems, being irregular and disky, could have 
high fractions of gas and dust. Simulations have 
shown that for high gas fractions, the outcome 
of a merger is sensitive to the treatment of the 
interstellar medium and could lead to the for- 
mation of rotatiqnally supported remnants (e.g., 



Robertson et al. 2006). Another scenario is that 



the strong star-formation episodes occur later in 
the merger timescales of z=2 ULIRGs than those 
of local ULIRGs due to differences in the dynam- 
ical configuration of their progenitors. 
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7. Conclusions 

We acquired Hubble Space Telescope NIC2 im- 
ages of 33 z ^2 ULIRGs to investigate the dynam- 
ical triggers of high-z, 24 /im-sclcctcd populations, 
revealed by Spitzer. Having a resolution clement 
of 1.3 kpc at 2=2, our ii-band images enable us 
to perform a comprehensive morphological analy- 
sis of these sources. Our findings are summarized 
as follows. 

1. The rest-frame optical luminosities of the 
sources in our sample are of the order 10 10 — 
10 11 L & and their mean half-light radius 
equals 2.66 kpc. On average, this observed 
R e g value is half of that of local ULIRGs, but 
surface brightness dimming might account 
for all or part of this difference. 

2. Taking into account binary systems, tidal 
tails and other irregular morphologies, the 
number of sources apparently involved in 
interactions is 17 out of 33. This is a 
lower limit since unresolved binaries, possi- 
ble faint companions, or diluted low surface- 
brightness features (of /z,y>22 mag arcsec~ 2 ) 
could exist in our sample. We simulated 
the effects of surface brightness dimming 
on our photometric analysis and found that 
extended structures such as tidal tails can 
typically be seen in systems as faint as 
n%H= 20. The mean apparent magnitude of 
the sources in our sample equals m# =20.89. 

3. Up to one- fifth of the sources in our sam- 
ple could be associated with minor mergers. 
The number of binary systems with compo- 
nents of comparable luminosity and signs of 
on-going interaction is only 2. This major- 
merger pair fraction is significantly smaller 
than that of the local ULIRG population. 

4. We performed a bulge-to-disk decomposition 
for all the sources that are brighter than 
niH—20. We found that these sources are 
disk-like with small or negligible bulges. 



always, suppressed in the compact systems. 
Its EW is high only in extended sources. The 
MIR obscuration correlates neither with the 
radial extent nor with the rest-frame optical 
luminosity of the sources. 

We conclude that a nonncgligiblc fraction of 
the sources in our sample do not resemble local 
ULIRGs, which are produced by major mergers 
and form elliptical galaxies. However, it is un- 
clear whether this result is characteristic of most 
z=2 ULIRGs or whether it relates to our selection 
criteria. In either case, it provides evidence that 
several interaction types can be associated with a 
ULIRG at z=2. 

K. D. wishes to thank Lee Armus and Lisa 
Storrie-Lombardi for useful discussions. This 
work is based on observations made with the 
NASA/ESA Hubble Space Telescope, obtained 
from the data archive at the Space Telescope Sci- 
ence Institute. STScI is operated by the Associ- 
ation of Universities for Research in Astronomy, 
Inc. under NASA contract NAS 5-26555. It is 
also based in part on observations made with the 
Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. Support 
for this work was provided by NASA through an 
award issued by JPL/Caltech. 



5. The AGN-dominated systems have the most 
compact and the PAH-dominated systems 
have the most extended morphology in the 
rest-frame optical images. The luminosity 
of the 7.7 jum PAH feature is often, but not 
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A. APPENDIX: Surface Brightness Dimming Simulations 

A question that naturally arises during the photometric analysis of intermediate and high z sources is 
what the effects are of surface br ightness dimming on the analysis performed. Although this issue has 
been addressed several times (e.g.. iGiavalisco et al. 19961 Hibbard h Vacca1ll997[ ). the results significantly 



depend on the type of sources under examination and the setup (e.g., central wavelength) of the observations. 
We therefore constructed a model to specifically address the following questions: (1) Up to which magnitude 
can secondary, faint nuclei be observed? (2) How do the observed magnitudes and effective radii relate to 
their intrinsic values? And (3) at which point does the bulgc-to-disk decomposition become unreliable? 

Our surface-brightness-dimming simulations use templates of local galaxies that are shifted to z=2. The 
templates were selected in a manner such that a large number of galaxy types are examined. We used archival 
data (Pis: A. S. Evans, R. Sharpies) of two ULIRGs, one LIRG, one quasar, and one elliptical galaxy in the 
local Universe observe d with ACS. Namely , these sources are IRAS 23128-5919, which is a binary ULIRG 
separated by 4.3 kpc ( Dasvra et al. I [2006ah . IC 4687, which is a LIRG with a significantly distorted disk, 



IRAS 09111-1007, which has secondary components, NGC 3156, and PG 1351+640. The sources were 
observed using either the 435W, the 475W, or the 625W filter, similar to the rest-frame wavelengths of our 
observations. 



To simulate z=2 galaxies from local Universe templates, we followed the prescriptions of IGiavalisco et al 



( 19961 ). First, we scaled down the angular size 9 of the template galaxy. The binning factor b was given by 

where Zi is the initial redshift of each template, Zf is the final redshift for all templates (which is equal to 
2), and s z is the pixel scale of the instrument used to perform the observations at a given redshift z. In 
our case, the instrumental pixel scales correspond to s Z/ =s(NICMOS)=0'/076 and Sz i =s(ACS)=0'/048. To 
ensure that the physical size of each source remains the same at all redshifts, its angular size ratio at Zf and 
Zi needs to be 

e H ~ {i+zi)*D L (z f ) ' 

where Dl(z) is the luminosity distance of a source at redshift z. 

We then accounted for the surface brightness dimming by multiplying the rescaled images with a dimming 
factor d that was equal to 

The factor a z is the constant that converts counts (or electrons) per second to flux. For the 160W filter of 
NICMOS, a Zf is 1.75 x 10~ 19 ergs cm" 2 A" 1 count- 1 . For ACS 435W, 475W, and 625W, the factor a Zi is 

equal to 3.14, 1.81, and 1.20 xl0~ 19 ergs cm -2 A 1 electron -1 , respectively. The quantity AX Z is the filter 
width for the instrument used for the observations at redshift z; it is equal to 4000, 1038, 1458, and 1442 
A for NICMOS 160W, ACS F435W, ACS, F475W, and ACS F625W, respectively. The factor x is used to 
match the intrinsic luminosity of the template with that of a source of a given magnitude at z—2. Given the 
observed magnitudes of our sources, we set the initial m^(init) value to either 19, 20, 21, or 22 mag. 

After convolving the resulting images with the NIC2 PSF, we added background noise to them. For this 
step, we computed and added to our images the median sky image of all data sets. To imitate the reduction 
of the real data, we then subtracted from these images the sky background of a random data set. We present 
the final simulated images in Figure [H We processed these images with sextractor and galfit in a manner 
similar to the NICMOS data. The magnitudes and effective radii, as they would be observed following the 
addition of the sky background, are denoted as m#(obs) and i? e ff(obs) and are presented in Table 5) 
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Already a visual inspection of the simulation results indicates that sources of m#(mit)=19 preserve to a 
large extent their initial light profiles. The sources of m_f/(init)=20 lose part of their low surface brightness 
components but still maintain an extended structure. Mergers of equal luminosity systems can be identified 
all the way up to m^r(init)=21 or an observed magnitude m#(obs) that is < 22.0. However, faint secondary 
nuclei are mostly seen up to m# (init)=20. For sources with m#(init) >21, one could often be unable to 
disentangle between minor mergers, single-nucleus remnants of major mergers, or quiescent sources. 

To quantify the surface-brightness-dimming effects on the photometric properties of our sources, we ran 
sextractor for all simulated galaxies. The results arc presented in Tabic 0] and depend on the initial light 
profile distribution. Compact distributions, such as PSFs and bulge components are more resilient to the 
surface brightness dimming than extended, disk-like distributions. For instance, we found that a spiral or 
irregular galaxy of TOff(init) =19 and 20 loses about 27% and 38% of its light, respectively, because of 
dimming effects. The observed effective radius is smaller by less than 14% and 25%, respectively, than its 
initial value. At m# (init) =21, about half of the light is lost, leading to R e g values that are underestimated 
by a factor of 2. Systems of mn (init) >22 are typically undetected at 3a levels. In contrast, the light losses 
are small for an elliptical galaxy. Its flux is decreased by 1% up to 22% and its radius is reduced by 2% up 
to 17% from its intrinsic values at m#(init) =19 and 22. 

We then aimed to address the reliability of the bulge-to-disk decomposition technique by examining 
how different its results are when shifting a local source to z=2, and when increasing (or decreasing) the 
brightness of a z—2 system by 1 mag. For this purpose, we ran galfit for some local sources and their 
simulated z=2 images, namely, the disk-dominated main component of IC 4687 and the elliptical galaxy 
NGC 3156. IRAS 23128-5919 and IRAS 09111-1007 are too perturbed to allow a decomposition. Compact 
light profiles such as those of PG 1351+640 were not analyzed either because they are degenerate between 
a PSF, and a small, unresolved bulge or disk. Since the light profiles of IC 4687 and NGC 3156 are close 
to those of a disk and a bulge, we examined the value of their Sersic index as a function of redshift and 
integrated magnitude. We find that the z=2 images of IC 4687 have n =1.05 and 1.06, respectively, for 
m#(init)=19 and 20. For comparison, the actual Sersic index of IC 4687 at a redshift of 0.0173 is 1.20. The 
z=2 images of NGC 3156 have n =5.67 and 5.28 for m H (init)=19 and 20, whereas at 0=0.004396, NGC 3156 
has n =6.70. The small differences for the elliptical galaxy are likely related to the spatial dilution of the 
initial (ACS) PSF with redshift. This comparison shows that the shape of the light profile can be accurately 
predicted at z—2. However, we wish to note that in these simulations, we used different magnitudes but 
identical initial structural parameters for the z=2, m#(init)=19 and 20 images of each source to properly 
examine the effects of surface brightness dimming. For observational data there is an extra difficulty; the 
best galfit solution can significantly depend on the initial parameter values. For the systems of m# (init)=19, 
a reasonable input file can be constructed by visual inspection of the images, e.g., by determining the size of 
the extended structures. For the systems of m#(init) >20, the surface brightness dimming does not always 
allow for the determination of an optimal input file. Galfit can then fall into a local (instead of the global) 
X 2 minimum. To avoid degenerate solutions, we opt to apply this decomposition technique only to sources 
that arc brighter than m#=20. 



11 



REFERENCES 

Alonso-Herrero, A., Rieke, G. H., Ricke, M, J., 
Colina, L., Perez- Gonzalez, P. G., & Ryder, S. 
D. 2006, ApJ, 650, 835 

Armus, L., Hcckman, T. M., & Miley, G. K. 1989, 
ApJ, 347, 727 

Armus, L, et al. 2007, ApJ, 656, 148 

Barker, E., ct al. 2006, NICMOS Instrument 
Handbook, Version 9.0 

Barnes, J. E. 1992, ApJ, 393, 484 

Bell, Eric F. et al. 2005, ApJ, 625, 23 

Bcrtin, E., & Arnouts, S. 1996, A&AS, 117, 393 

Brandl, B. R. ct al. 2006, ApJ, 653, 1129 

Bridge, C. R., et al. 2007, ApJ, 659, 931 

Bushousc, H. A., ct al. 2002, ApJS, 138, 1 

Chapman, S. C, Smail, I., Windhorst, R., 
Muxlow, T., & Ivison, R. J. 2004, ApJ, 611, 
732 

Conselice, C. J. 2006, ApJ, 638, 686 

Dasyra, K. M., ct al. 2006a, ApJ, 638, 745 

Dasyra, K. M., et al. 2006b, ApJ, 651, 835 

Dickinson, M., Sosey, M., Rieke, M., Bohlin, R. 
C, & Calzetti, D. 2002, hstc, conf, 231 

Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 

Downcs, D., & Solomon, P.M. 1998, ApJ, 507, 615 

Dressier, A., Lyndcn-Bcll, D., Burstein, D., 
Davies, R. L., Faber, S. M., Terlevich, R., & 
Wegner, G. 1987, ApJ, 313, 42 

Elbaz, D., Cesarsky, C. J., Chanial, P., Aussel, 
H., Franceschini, A., Fadda, D., & Chary, R. 
R. 2002, A&A, 384, 848 

Fadda, D., et al. 2006, AJ, 131, 2859 

Fazio, G., et al. 2004, ApJS, 154, 10 

Flores, H. et al. 1999, ApJ, 517, 148 

Genzel, R., et al. 1998, ApJ, 498, 579 



Genzel, R., Tacconi, L. J., Rigopoulou, D., Lutz, 
D., & Tecza, M. 2001, ApJ, 563, 527 

Giavalisco, M., Livio, M., Bohlin, R. C, Mac- 
chetto, F. D., Stccher, T. P. 1996, AJ, 112, 369 

Hammer, F., Flores, H., Elbaz, D., Zheng, X. Z., 
Liang, Y. C, & Cesarsky, C. 2005, A&A, 430, 
115 

Haussler, B. ct al. 2007, ApJS, 172, 615 

Hibbard, J. E., & Vacca, W. D. 1997, AJ, 114, 
1741 

Houck, J., et al. 2004, ApJSS, 154, 18 

Ishida, C. M. 2004, Ph. D. Thesis, University of 
Hawaii 

Joseph, R. D., & Wright, G. S. 1985, MNRAS, 
214, 87 

Kartaltepe, J. S., et al. 2007, ApJ, submitted 
(arXiv/0705.2266vl) 

Kim, D.-C, Veilleux, S., & Sanders, D. B. 2002, 
ApJS, 143, 277 

Kormendy, J., Sanders, D. B. 1992, ApJ, 390L, 53 

Krist, J. E., & Hook, R. N. 1997, HST Calibration 
Workshop, Space Telescope Science Institute, 
S. Casertano eds. 

Lacy, M., et al. 2005 ApJS, 161, 41 

Larson, R. B., & Tinsley, B. M. 1978, ApJ, 219, 
46 

Le Floc'h, E. et al. 2005, ApJ, 632, 169 

Melbourne, J., Koo, D. C, & Lc Floc'h, E. 2005, 
ApJ, 632, 65L 

Metcalfe, N., Shanks, T., Weilbacher, P. M., Mc- 
Cracken, H. J., Fong, R., & Thompson, D. 2006, 
MNRAS, 370, 1257 

Murphy, T. W., Jr., Soifer, B. T., Matthews, K., 
& Armus, L. 2001, ApJ, 559, 201 

Murphy, T. W., Jr., Armus, L., Matthews, K., 
Soifer, B. T., Mazzarella, J. M., Shupe, D. L., 
Strauss, M. A., & Neugcbauer, G. 1996 AJ, 111, 
1025 



12 



Peng, C. Y., Ho, L. C, Impey, C. D., & Rix, H.-W. 
2002, AJ, 124, 266 

Perez-Gonzalez, P. G., et al. 2005, ApJ, 630, 82 

Rigopoulou, D., Spoon, H. W. W., Genzel, R., 
Lutz, D., Moorwood, A. F. M., & Tran, Q. D. 
1999, AJ, 118, 2625 

Robertson, B., Bullock, J. S., Cox, T. J., Di Mat- 
teo, T., Hernquist, L., Springel, V., & Yoshida, 
N. 2006 ApJ, 645, 986 

Sajina, A., Yan, L., Armus, L., Choi, P., Fadda, 
D., Helou, G., & Spoon, H. 2007, ApJ, 664, 713 

Sajina, A., et al. 2008, ApJ, submitted 

Sanders, D. B., Soifer, B. T., Elias, J. H., Madore, 
B. F., Matthews, K., Neugebauer, G., & Scov- 
ille, N. Z. 1988a, ApJ, 325, 74 

Sanders, D. B., Soifer, B. T., Elias, J. H., Neuge- 
bauer, G., & Matthews, K. 1988b, ApJ, 328, 
35L 

Scoville, N. Z., et al. 1998, ApJ, 492, L107 

Scoville, N. Z., et al. 2000, AJ, 119, 991 

Simmons, B. P., fc Urry, C. M. 2008, ApJ, sub- 
mitted, (|arXiv:0804. 1363vl) 

Soifer, B. T., Neugebauer, G., Matthews, K., 
Egami, E., & Weinberger, A. J. 2002, AJ, 124, 
2980 

Spoon, H. W. W., Marshall, J. A., Houck, J. R., 
Elitzur, M., Hao, L., Armus, L., Brandl, B. R., 
& Charmandaris, V. 2007, ApJ, 654, 49L 

Springel, V., Di Matteo, T., & Hernquist, L. 2005 
ApJ, 620, L79 

Springel, V., & Hernquist, L. 2005 ApJ, 622, L9 

Surace, Jason A., Sanders, D. B., Vacca, W. D., 
Veilleux, S, & Mazzarella, J. M. 1998, ApJ, 492, 
116 

Surace, J. A., & Sanders, D. B. 2000, AJ, 120, 604 

Surace, J. A., Sanders, D. B., & Evans, A. S. 2000, 
ApJ, 529, 170 

Swinbank, A. M., Chapman, S. C, Smail, I., Lind- 
ner, C, Borys, C, Blain, A. W., Ivison, R. J., 
& Lewis, G. F. 2006, MNRAS, 371, 465 



Tacconi, L. J., Genzel, R., Lutz, D., Rigopoulou, 
D., Baker, A. J., Iserlohe, C, & Tecza, M. 2002, 
ApJ, 580, 73 

Tacconi, L. J. et al. 2006, ApJ, 640, 228 

Taylor, V. A., Jansen, R. A., Windhorst, R. A., 
Odewahn, S. C, & Hibbard, J. E. 2005, ApJ, 
630, 784 

Tran, Q. D., et al. 2001, ApJ, 552, 527 

Veilleux, S., Kim, D.-C, Sanders, D. B., Maz- 
zarella, J. M., & Soifer, B. T. 1995, ApJS, 98, 
171 

Veilleux, S., Kim, D.-C, & Sanders, D. B. 1999, 
ApJ, 522, 113 

Veilleux, S., Kim, D.-C, & Sanders, D. B. 2002, 
ApJS, 143, 315 

Veilleux, S., et al. 2006, ApJ, 643, 707 

Yan, L., et al. 2004, ApJS, 154, 75 

Yan, L., et al. 2005, ApJ, 628, 604 

Yan, L., et al. 2007, ApJ, 658, 778 

Zheng, X. Z., Hammer, F., Flores, H., Assemat, 
F., & Pelat, D. 2004, A&A, 421, 847 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 



13 



Table 1 

Source List and Observation Log 





RA a 


De 








e b 


log(Z-i 4 / Lq) c log( .L ir /£q) 




MIR b 


t d 
tint 






(2000) 


(2000) 
















type 


(s) 


(A) 


MIPS42 


17:17 


58.54 


+ 59:28 


16 


18 


1 


95 


12 


68 








5120 


4750-6100 


MIPS78 


17:15 


38.18 


+ 59:25 


40 


08 


2 


65 


12 


76 






mixed 


5248 


3840-4930 


MIPS180 


17:15 


43.73 


+ 58:35 


32 


97 


2 


47 


12 


57 






obscured 


5248 


4040-5190 


MIPS227 


17:14 


55.99 


+ 58:38 


16 


44 


1 


87 


12 


39 






mixed 


2560 


4880-6270 


MIPS289 


17:13 


50.04 


+ 58:56 


54 


45 


1 


86 


11 


80 


12.9 




PAH 


2560 


4900-6300 


MIPS464 


17:14 


39.60 


+ 58:56 


32 


07 


1 


85 


12 


16 






mixed 


5120 


4910-6320 


MIPS506 


17:11 


38.52 


+ 58:38 


38 


58 


2 


52 


12 


74 






mixed 


2560 


3980-5110 


MIPS8196 


17:15 


10.17 


+ 60:09 


54 


54 


2 


60 


12 


53 






obscured 


2688 


3890-5000 


MIPS8242 


17:14 


33.19 


+ 59:39 


11 


31 


2 


45 


12 


74 






mixed 


2560 


4060-5220 


MIPS8245 


17:15 


36.31 


+ 59:36 


14 


73 


2 


70 


12 


55 






mixed 


5248 


3780-4870 


MIPS8327 


17:15 


35.62 


+ 60:28 


24 


49 


2 


48 


12 


40 






mixed 


5248 


4020-5170 


MIPS8342 


17:14 


11.45 


+60:11 


09 


20 


1 


57 


11 


93 


12.7 




AGN 


2688 


5450-7000 


MIPS8493 


17:18 


04.99 


+ 60:08 


32 


28 


1 


80 


11 


49 






PAH 


2688 


5000-6430 


MIPS15840 


17:19 


22.44 


+60:05 


00 


42 


2 


30 


12 


50 






AGN 


5248 


4240-5450 


MIPS15880 


17:21 


19.61 


+ 59:58 


17 


87 


1 


64 


12 


51 


12.9 




obscured 


2560 


5310-6820 


MIPS15928 


17:19 


17.47 


+60:15 


19 


62 


1 


52 


12 


04 






mixed 


2688 


5560-7140 


MIPS15949 


17:21 


09.21 


+ 60:15 


01 


66 


2 


15 


12 


45 






AGN 


2688 


4440-5710 


MIPS15958 


17:23 


24.21 


+ 59:24 


55 


62 


1 


97 


12 


28 






mixed 


2560 


4710-6060 


MIPS15977 


17:18 


55.66 


+ 59:45 


45 


61 


1 


85 


12 


28 


13.0 




AGN 


2560 


4910-6320 


MIPS16059 


17:24 


28.34 


+60:15 


33 


08 


2 


43 


12 


46 






mixed 


5376 


4080-5250 


MIPS16080 


17:18 


44.83 


+ 60:01 


15 


82 


2 


04 


12 


18 






obscured 


2688 


4610-5920 


MIPS16095 


17:23 


59.77 


+ 59:57 


52 


42 


1 


81 


12 


12 






mixed 


5120 


4980-6410 


MIPS16113 


17:21 


26.47 


+60:16 


46 


27 


1 


90 


11 


88 






obscured 


2688 


4830-6210 


MIPS16122 


17:20 


51.50 


+60:01 


48 


75 


1 


97 


12 


13 






mixed 


5376 


4710-6060 


MIPS16144 


17:24 


22.03 


+ 59:31 


50 


63 


2 


13 


12 


01 






PAH 


2560 


4470-5750 


MIPS22204 


17:18 


44.40 


+ 59:20 


00 


96 


2 


08 


12 


63 






obscured 


2560 


4550-5850 


MIPS22277 


17:18 


26.73 


+ 58:42 


42 


12 


1 


77 


12 


36 






mixed 


2560 


5050-6500 


MIPS22303 


17:18 


48.87 


+ 58:51 


14 


87 


2 


34 


12 


57 






obscured 


5248 


4190-5390 


MIPS22530 


17:23 


03.31 


+ 59:16 


00 


55 


1 


96 


12 


10 


13.0 




PAH 


2560 


4730-6080 


MIPS22558 


17:20 


45.15 


+ 58:52 


21 


50 


3 


20 


12 


88 






AGN 


5248 


3330-4290 


MIPS22651 


17:19 


26.54 


+ 59:09 


29 


12 


1 


73 


11 


86 


12.9 




mixed 


2560 


5130-6590 


MIPS22661 


17:18 


19.63 


+ 59:02 


42 


68 


1 


75 


11 


93 






AGN 


5248 


5090-6550 


MIPS22699 


17:20 


47.52 


+ 59:08 


14 


71 


2 


59 


12 


17 






mixed 


5248 


3890-5010 


a Units of right ascon 




hours, n 




id 




nds, and 


units of declin 


ation arc degre 




, and ar 


cseconds. The 


coordinates of the sour 




those of the 


IRAC 3.6 fj.i 


n catalog 


of Lacy ct al 


2005). 










b The redshift 


s and the MIR 




sed for the type 


classifications arc fron 


Ibaiina et a 


• 


1200711. The spectral type nomencla- 


turc adopted he 


rc simply dene 


tes the d 








npor 


cnt of th 


c MIR spectn 


m, namely the 


AGN continu 


am, the 


PAH emission 


features , or the 


SiO absorptio 






9.7 fj. 




If the domin 


ant componei 


t is not obv 




s, we simply describe 


it as a mixed 
































c Thc IR and 


14 {Am 


uminos 




n Saji 




ct a 


. (2008) 


and Sajina et 


al. | 12001) r 


espectively. The 


14 ^m 




be considered a 


lower 


imit of 


L IR whe 


n the latter is n 


ot availa 


ble. 












d Total integration t 


me of the NIC2 


3xp 






















e Rest-frame a 


vavelen 


gth ran 


of the 


obscrvat 




s pc 


--formed with the NIC2 




th 


e F160W filto 
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Table 2 

Global Photometric Parameters 



Galaxy 








^cff a vL u (opt) a interaction 








(mag arcscc ) 


(kpc) ( L q) signs 



MIPS42 


21 


79 


-24. 


10 


21 


29 


2 


66 


2 


68X10 10 




MIPS78 


21 


.98 


-24 


72 


21 


95 


3 


13 


4 


72X10 10 


none'* 


MIPS180 


22 


.15 


-24 


37 


20 


85 


1 


77 


3 


41 X 10 10 


noncb 


MIPS227 


18 


.85 


-26 


93 


18 


32 


2 


63 


3 


62X10 11 




MIPS289 


19 


72 


-26 


04 


20 


01 


3 


84 


1 


60X10 11 


distorted 


MIPS464 


21 


45 


-24. 


30 


20 


62 


2 


29 


3 


19X 10 10 


distorted* 1 


MIPS506 


21 


.66 


-24. 


91 


20 


60 


1 


.97 


5 


60x 10 10 


distorted 


MIPS8196 


18 


.68 


-27 


97 


18 


76 


3 


31 


9 


.46X 10 11 


distorted 


MIPS8242 


20 


21 


-26 


28 


20 


.85 


4 


34 


1 


.99X 10 11 


distorted* 3 


MIPS8245 


22 


.21 


-24. 


54 


21 


.17 


1 


96 


3 


98X 10 10 




MIPS8327 


21 


06 


-25. 


47 


20 


09 


2 


.07 


9 


35X 10 10 


distorted 


MIPS8342 


21 


.54 


-23 


77 


20 


68 


2 


28 


1 


96X 10 10 


none*- 1 


MIPS8493 


21 


23 


-24. 


45 


21 


18 


3 


29 


3 


.68X 10 10 


distorted 


MIPS15840 


21 


67 


-24. 


66 


20 


96 


2 


36 


4 


.44X 10 10 


none d 


MIPS15880 


21 


26 


-24 


17 


21 


74 


4 


22 


2 


83X 10 10 


distorted* 1 


MIPS15928 


19 


65 


-25. 


57 


19 


31 


2 


.88 


1 


.04X 10 11 


distorted 


MIPS15949 


20 


.73 


-25 


42 


20 


00 


2 


36 


9 


00X 10 10 


none 


MIPS15958 


21 


09 


-24. 


83 


19 


25 


1 


43 


5 


19X 10 10 


none 


MIPS15977 


19 


83 


-25. 


92 


19 


05 


2 


35 


1 


.42X 10 11 


distorted 


MIPS16059 


20 


49 


-25. 


98 


20 


.19 


2 


82 


1 


51 X 10 11 


binary 


MIPS16080 


20 


31 


-25. 


70 


19 


62 


2 


43 


1 


16X 10 11 


distorted 


MIPS16095 


19 


.96 


-25. 


73 


18 


93 


2 


09 


1 


20X 10 11 




MIPS16113 


21 


40 


-24 


42 


20 


94 


2 


71 


3 


.58X10 10 


distorted 


MIPS16122 


21 


.36 


-24 


56 


21 


.14 


3 


03 


4 


08X 10 10 


distorted 


MIPS16144 


20 


.85 


-25. 


27 


20 


80 


3 


24 


7 


83X 10 10 


distorted 6 


MIPS22204 


20 


17 


-25. 


89 


18 


50 


1 


.53 


1 


39X 10 11 


none 


MIPS22277 


20 


.08 


-25. 


55 


19 


.10 


2 


15 


1 


Olx 10 11 


none b 


MIPS22303 


22 


.76 


-23 


61 


22 


.13 


2 


44 


1 


70X 10 10 




MIPS22530 


21 


.07 


-24 


83 


21 


88 


1 


.87 


5 


25X 10 10 


binary 


MIPS22558 


21 


62 


-25. 


58 


21 


39 


2 


71 


1 


03X 10 11 


distorted 


MIPS22651 


20 


11 


-25. 


46 


19 


32 


2 


35 


9 


35X 10 10 


none* 1 


MIPS22661 


20 


22 


-25 


38 


19 


48 


2 


40 


8 


65 X 10 10 


none 


MIPS22699 


22 


.19 


-24. 


45 


20 


88 


1 


75 


3 


.67X 10 10 


none 



a The /f-band (or rcstframe optical) apparent and absolute magnitudes, half-light radii, 
and surface brightnesses were derived using Scxtractor, which models the galaxies as series of 
isophotal ellipses. For the spatially overlapping binary systems, these numbers arc computed 
for the system as a whole. No spectroscopic fc-corrcctions have been applied to these quantities. 

^In these systems the nucleus docs not seem distorted, but secondary low surface-brightness 
sources exist within a projected separation of 10 kpc. Although these systems may be dy- 
namically related, no prominent indication of interaction between them can be seen at our 
surface- brightness limits. 

c In these cases, the bright nucleus has a faint component in the same envelope that could 
cither be a secondary nucleus or a bright star- formation knot. 

d Bright sources of Am <1.2 exist close to these sources but there are no signs of interaction 
between them. 

e Secondary, low surface- brightness structures exist around the distorted nucleus of these 
systems. 
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Table 3 

Bulge-to-disk Decomposition Parameters 



Galaxy 


m H (bulge) a 


B eff (bulge) B 
(kpc) 


m H (disk) a 


B eff (disk) a 
(kpc) 


bulgc-to-disk 


rcsiduals b 

% 


c 




Reff(n) 
(kpc) 


m^(PSF) <= 


MIPS227 


19.90 


2.37 


19.30 


1.98 


0.57 


6.25 


1.23 


18.83 


2.35 


21.97 


MIPS289 


20.89 


1.25 


20.14 


5.76 


0.50 


3.51 


0.60 


20.07 


3.66 


21.91 


MIPS8196 


19.15 


2.06 


19.64 


4.36 


1.57 


0.97 


2.65 


18.64 


3.10 


22.14 


MIPS15928 


22.60 


1.20 


19.78 


2.64 


0.07 


14.2 


0.67 


19.85 


2.57 


21.49 


MIPS15977 


20.86 


1.13 


20.21 


2.01 


0.55 


1.00 


0.54 


20.09 


2.08 


21.31 


MIPS16095 


21.15 


1.38 


20.36 


2.03 


0.48 


0.72 


1.35 


20.41 


1.94 


23.31 



a The best-fit solutions for the bulge and the disk magnitudes and half-light radii were derived using galfit. For an exponential disk, the half-light radius 
is related to the scale length as R c ff= 1-68-R^. 

'-'The residuals arc given as a percentage fraction of the initial flux. 

c For comparison with the bulge-to-disk decomposition results, we also perform a free Sersic index n versus PSF decomposition. We provide the value of 
n and the magnitude of the Sersic component, as well as the PSF magnitude. 
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Table 4 

Photometric parameters of simulated z=2 galaxies 



Galaxy 




init) a 


< ^ H (init) > 


R eff (init) a 




obs) b 


< Mff (obs) > 


*eff <° b3 > b 


template 






(mag arcscc 2 ) 


(kpc) 






(mag arcscc - 2 ) 


(kpc) 


IC 4687 


19 


00 


20.61 


7.01 


19 


37 


20.65 


6.03 


IC 4687 


20 


.00 


21.61 


7.01 


20 


.51 


21.49 


5.24 


IC 4687 


21 


.00 


22.62 


7.01 


21 


84 


21.90 


3.43 


IC 4687 c 


22 


.00 


23.62 


7.01 










IRAS 09111-1007 


19 


.00 


20.37 


6.29 


19 


31 


20.45 


5.65 


IRAS 09111-1007 


20 


.00 


21.37 


6.29 


20 


.53 


21.30 


4.77 


IRAS 09111-1007 


21 


00 


22.37 


6.29 


21 


74 


22.13 


4.00 


IRAS 09111-1007 c 


22 


.00 


23.37 


6.29 










IRAS 23128-5919 


19 


.00 


19.70 


4.62 


19 


20 


19.63 


4.06 


IRAS 23128-5919 


20 


.00 


20.70 


4.62 


20 


.44 


20.36 


3.22 


IRAS 23128-5919 


21 


.00 


21.70 


4.62 


21 


.60 


21.30 


2.91 


IRAS 23128-5919 


22 


.00 


22.70 


4.62 










NGC 3156 


19 


.00 


17.88 


1.99 


19 


.06 


17.82 


1.89 


NGC 3156 


20 


.00 


18.88 


1.99 


20 


10 


18.78 


1.82 


NGC 3156 


21 


.00 


19.88 


1.99 


21 


.17 


19.74 


1.73 


NGC 3156 


22 


.00 


20.88 


1.99 


22 


26 


20.74 


1.66 


PG 1351+640 


19 


.00 


17.65 


1.79 


19 


.12 


17.62 


1.67 


PG 1351+640 


20 


00 


18.65 


1.79 


20 


17 


18.56 


1.59 


PG 1351+640 


21 


.00 


19.65 


1.79 


21 


.21 


19.50 


1.52 


PG 1351 + 640 


22 


.00 


20.65 


1.79 


22 


25 


20.45 


1.50 



a Magnitudc and effective radius of each redshifted source prior to the addition of the sky background. 

b Magnitude and effective radius of each redshifted source as measured after the addition of an appropriate sky back- 
ground. All quantities arc measured with sextractor, as for the NICMOS observations. 

c These sources cannot be detected at 3tr levels. 
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8 




Fig. 1. — Rcdshift distribution of the ULIRG sample observed with the Hubble Space Telescope NIC2 camera. 




Fig. 2. — Reduced NIC2 images of the z~2 ULIRGs in our sample. All the images in this figure correspond 
to a 3'.' 8 x 3'.' 8 field of view and are displayed in a squared scale. 
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Fig. [2]- continued. 
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Fig. [2]- continued. 
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Fig. 3.— Histograms of the apparent magnitude distribution {top left), the effective radius distribution {top 
right), the R — H colors {bottom left), and the H— IRAC(3.6 /an) colors {bottom right) of the sources in our 
sample. 
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Fig. 4. — Bulge-to-disk decomposition results for the bright, extended sources in our sample. For each source, 
the first panel corresponds to the original image viewed in a 1'.'9 x 1'.'9 frame. The second panel corresponds 
to the residuals after the subtraction of the disk component. The final, bulge-subtracted residuals are 
presented in the third column. Whereas the decomposition steps are presented here individually for the sake 
of clarity, the actual galfit best-fit solution was found by simultaneously solving for the optimal bulge and 
disk parameters. 
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Fig. 5. — Observed ff-band (rest-frame optical) photometric properties as a function of the MIR spectral 
type. The apparent magnitudes (left) and the luminosities (right) are plotted vs. the effective radii for AGNs 
(stars), starbursts (boxes) and obscured systems (crosses). In this diagram, we are not using sources with 
composite IRS spectra, i.e., sources with simultaneous PAH emission and SiO absorption features, and AGN 
continua (see Table [lj. 
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Fig. 6. — Luminosity (left) and equivalent width (right) of the 7.7 fxm PAH feature as a function of the 
effective radius. All sources with 7.7 fim feature detections are used in this diagram, regardless of the 
dominant component of the MIR spectrum. Filled circles correspond to sources with detections and open 
circles correspond to sources with upp er limits on their 7.7 /zm luminosities. The dashed line on the right 
panel corresponds to the threshold that lSaiina et al." ( 2007 ) set to discriminate between sources with strong 
and weak PAH emission. 
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Fig. 7. — MIR optical depth, as measured from the 9.7 /im SiO feature, is plotted as a function of the rest- 
frame optical effective radius (left) and the rest-frame optical luminosity (right). Filled circles are sources 
with detections and open circles are sources with lower limits on their MIR obscuration. 
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Fig. 8. — Results from a surface-brightness dimming model that simulates how local galaxies of different 
morphological types would be observed at z—2. The first column shows the 40" x 40" ACS images of the 
local template galaxies. The next four columns show the simulated NIC2 images of the same galaxies that 
are scaled to match the intrinsic luminosity of a z=2 source with an apparent magnitude of mjj(init)=19, 
20, 21, and 22 respectively, prior to its dilution into the sky background. The magnitudes of these objects 
as observed after the addition of the background are tabulated in Table [U The size of the NIC2 simulated 
images corresponds to 3"8 x 3"8, similar to the observations in Fig. [2] 
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